{
j++
for ( i=1 ; i<=NF ; i++ ){
    av[j] += $i
    pt[i,j] = $i
    }
av[j] /= NF
nf = NF
}

END {
   for ( j=1 ; j<=NR ; j++ ) {
      sig[j] = 0.0 
      for ( i=1 ; i<=nf ; i++ ) {
         sig[j] += (pt[i,j]-av[j]) * (pt[i,j]-av[j])
       }
   sig[j] = sqrt(sig[j]/(NR-1))
   print av[j],sig[j]
   }
}
